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We study T — spin density wave transitions in two dimensional Fermi liquids in which the 
ordering wavevector Q is such that the tangents to the Fermi line at the points connected by Q are 
parallel (e.g. Q — 2pF in a system with a circular Fermi line) and the Fermi line is not flat. We 
show that the transition is first order if the ordering wave vector Q is not commensurate with a 
reciprocal lattice vector, G, i.e. Q 7^ G/2. If Q is close to G/2 the transition is weakly first order 
and an intermediate scaling regime exists; in this regime the 2pF susceptibility and observables such 
as the NMR rates Ti and T2 have scaling forms which we determine. 
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I. INTRODUCTION 

Quantum phase transitions have attracted substantial 
recent interest. Aatiferromagnet-singlet transitions in 
insulating magnetsEra and ferromagujEitic and antiferro- 
magnetic transitions in Fermi liquidsB'El have been stud- 
ied in detail, and the crossover between the insulating 
and fermi liquid critical points in two spatial dimensions 
has also been studiedQ. Here we consider an important 
case which has not so far been discussed in the literature, 
namely what we call the "2pf^" spin or charge density 
wave transition of a fermion system. By ^^2pp^^ we mean 
an ordering wavevector Q which connects two paints on 
the Fermi line with parallel tangents [See Fig. 1]eI. For a 
circular Fermi line any vector Q of magnitude 2pp con- 
nects two such points. In this paper we consider explicitly 
the spin density wave case, but our results can be applied 
with only minor modifications to the charge density wave 
case. We assume that the Fermi line is not straight. We 
also assume that Fermi liquid theory adequately repre- 
sents the non-critical properties of the fermions. If it does 
not, our results do not apply. We briefly discuss one non- 
Fermi liquid scenario in the conclusion. One important 
motivation for studying the 2pF case is the high Tc su- 
perconducting material La2-xSrxCuOi, iruwhich strong 
magnetic fluctuations have been observecH; the fluctu- 
ations are peaked at an x-dependent wavevector Q{x) 
which is claimed to be a "2pi?" wavevector of the Fermi 
line calculated, by standard band-structure techniques for 
this materiaLj. Our results may also be relevantJpr quasi 
two dimensional materials such as TTF-TCNQU. 

In order to study critical phenomena analytically, one 
expands about a mean field solution. If Fermi liquid the- 
ory is a good starting point, then the appropriate mean 
field theory is the Random Phase Approximation (RPA), 



in which the susceptibility q) is given in terms of the 
interaction constant g and the polarizability of noninter- 
acting fermions, no(w,q) by x(c^,q) = no(w,q)/(l - 
17^110(0;, q)). The transition occurs at the wavevector Q 
as g is increased to the point that (7^X10(0, Q) = 1. In 
three spatial dimensions, no(0, q) is not maximal at any 
2pp wavevector, moreover dIV[j{Q, q) / dq is logarithmically 
divergent as q 2pp. Therefore a "2pp" transition is 
impossible in c? = 3 and so we focus on d = 2 in this 
paper. 

In two dimensional case, no(0, q) is so strongly peaked 
at q = 2pp that it is natural to assume that the spin 
density instability happens at "2^^"- The ^^2pp" case 
is difficult to treat by the methods used Dixviously to 
study phase transitions at other momentaQ'El. In these 
works the fermions are "integrated out" and the prob- 
lem is reduced to a model of interacting bosonic spin 
fluctuations. In the "2pp" case the action functional ob- 
tained by integrating out the fermions has coefficients 
which are singular and non-analytic because the fermion 
response functions are non- analytic for Q — 2pp. These 
non-analyticities lead to divergences in the action as 
T ^ ajid make it difficult to apply the conventional 
approachoQ. Instead, in this paper we apply a perturba- 
tive renormalization group technique to a model which 
includes both spin fluctuations and fermions. 

Our perturbation parameter is 1/iV, the fermion spin 
degeneracy. The leading order of the perturbation theory 
is the familiar RPA approximation. The next order is a 
theory of electrons interacting by exchanging RPA fluc- 
tuations. We show that this theory is infrared divergent. 
We sum the leading infrared divergent contributions us- 
ing the renormalization group. 

The behavior of spin fluctuations changes dramati- 
cally if their wave vectors are close to half of a recip- 
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rocal lattice vector, G. The important parameter is 
AG = IQ — G/2|. If AG is sufficiently small we must dis- 
tinguish two regimes in the renormalization group flow: 
large momenta, where the infrared cutoff is greater than 
AG and small momenta, where it is less. For large mo- 
menta the divergences are logarithmic; the logarithms 
may be summed by renormalization group to power laws 
and we use the 1/iV-expansion to find the exponents. 
Although the physical value of = 2, the small value 
of the numerical coefficients in front of these logarithms 
suggests that the exponents obtained in the first order 
in 1/N are close to their exact values at iV = 2. For 
small momenta, the divergences are much stronger and, 
we show, drive the transition first order as soon as the 
regime of small momenta is reached. 




FIG 1. Sketch of Fermi line and important wavevec- 
tors. The Fermi line shown here is similar to that claimed 
to be appropriate to Lai.^Sr ij^CuOA- The ordering 
wavevector Q for spin fluctuations connects two points 
on the Fermi line. It is assumed that the tangent to the 
Fermi line at one end of the vector Q is parallel to the 
tangent to the Fermi line at the other end. We have also 
shown a typical momentum of a spin fluctuation q. We 
parametrize the vector k = q — Q by its cartesian com- 
ponents and k\\ in the coordinate system (shown by 
dashed lines) associated with the Fermi line at the points 
connected by Q. 



The outline of this paper is as follows. In Section II we 
deflne the model, derive the RPA theory and make a con- 
venient scaling of variables. In Section III we analyze the 
fluctuation corrections and derive renormalization group 
equations in the regime of large momenta. In Section 
IV we derive analogous equations in the regime of small 



momenta and show that they imply that the transition 
is first order. In Section V we discuss the physical conse- 
quences of our results. Section VI contains a summary of 
the results, a discussion of their relation to previous work 
on quantum critical phenomena and correlated electrons, 
a note on the extension of our results to a charge density 
wave transition and a conclusion. 



II. MODEL AND RANDOM PHASE 
APPROXIMATION 

Our starting point is a Hamiltonian, i7, describing 
fermions moving in a lattice and interacting with each 
other via a short range four fcrmion interaction W: 

pa p.p' ,q.ct,(3 

(1) 

We assume that a T = spin density wave transition 
to a state with long range order at wave vector Q occurs 
as W is increased to a critical value Wc- Because we ex- 
pect the physics in this region to be determined by the 
exchange of spin density fluctuations we use a Hubbard- 
Stratonovich transformation to recast Eq. (|l|) as a theory 
of fermions coupled to spin fluctuations Sq. The theory 
is described by the action 

pa 

+ J2D-\u;,q)S^,gS.^,g (2) 

p.q.e,uj 

Here G{e,p) is the fermion Green function, D{uj,q) is 
the spin fluctuation propagator and gi, is a bare coupling 
constant derived from W. When Hubbard-Stratanovich 
transformation is applied to Eq. (|l|), the result is 
Di,{uj,q) = 1, = W and Gb{€,p) is the non-interacting 
fermion Green function, i.e. 



The interaction between spin fluctuations and fermions 
changes the form of the fermion Green function and spin 
fluctuation propagator. We assume that the effects of the 
short scale fluctuations which do not become singular at 
the critical point can be described by conventional Fermi 
liquid renormalizations. 

The leading order of perturbation theory for the action 
is the Random Phase Approximation (RPA) which 
takes into account the renormalization of the spin prop- 
agator by the electron polarization bubble, DQ^{uj,q) = 
D^\u;,q)-UoiiJ,q). 



2 



We shall be interested in momenta close to the mo- 
mentum Q at which no(0,g) is maximal. For wavevec- 
tors near Q the momentum and frequency dependences 
of no(0, q) are non-analytic and controlled by Fermi line 
singularities. Because the singular behavior of no(w,(j) 
is controlled by the distance from q to the Fermi line, 
it will be convenient to parametrize the momentum q in 
terms of the variables fcy and k±_ shown in Fig. 1. 

The fermion polarizability no(w, g) can be calculated 
by summing all diagrams which are irreducible with re- 
spect to the fermion-fermion interaction and have two 
external ^i^^q legs. This generalizes the RPA by includ- 
ing Fermi liquid corrections. This sum has contributions 
from short length scale processes which give Ho (07,(7) 
an analytic dependence on q and uj and also contribu- 
tions from Fermi line singularities, which lead to a non- 
analytic dependence of no(w,(7) on q and lo. Thus we 
write no(w,q) = n"'°'(tj,g) + n^™9(w,q), within Fermi 
liquid theory singularities come from the diagram shown 
in Fig. 2. The analytic expression corresponding to this 
diagram is 



G{e + io,p + q)G{e,p) (4) 



where 50 represents the interaction constant renormal- 
ized by Fermi liquid corrections. To obtain the form of 
Fermi line singularities we expand the spectrum of the 
fermions in the vicinity of the points ±Q/2, obtaining 
for the fermion Green function 



G(e,p) 



1 



'JFP\\ 



2po 



(5) 



Here Wi? is the renormalizcd Fermi velocity, po is the ra- 
dius of curvature of the Fermi line and p\\ {p±) are momen- 
tum components normal (tangential) to the Fermi line as 
measured from the points ±Q/2 and z is the quasiparticle 
residue. Note that coordinates p\\ and p±^ are compatible 
with the spin fluctuation coordinates fc||, fcj,. 



All dominant infrared contributions come from pro- 
cesses in which an electron is scattered from one small re- 
gion of fermion momenta around Q/2 to another around 
— Q/2. It is convenient to introduce dimensionless mo- 
menta parameterizing these regions and rescale all fields 
so that the resulting action does not contain dimensionful 
variables. We choose 



Pi- ^/PqPfPa 

P\\ PFP\\ 

e vfPf£ 



Suj.k 

90 



-1 -5/2 -1/2 
'^F Pf Pq ^<^,p 

-1 -1/2 -3/2 r, 
5o Pf ^<^.k 



(6) 



After this transformation we may assume we are dealing 
with a circular Fermi surface of unit radius but with an 
upper cut off on angular integrals of the order of pf /po 
and an upper cut off on radial integrals of the order of 1; 
the dimensionless parameter controlling the incommen- 
surability becomes 



AG 



M 

2pF 



In the new variables the action retains the general form 
a) but the Green function of the fermions changes to 



G{e,p) 



1 



ize-pii +p\/2 
while the bare spin fluctuation propagator becomes 



(7) 



Db{uj,q) = g 



VF 



Evaluating the diagram shown in Fig. 2 yields the sin- 
gular part of the fermion polarizability in d = 2 



n^''™^ ) (o;, q) = i?e 
27rz 



-11 



(8) 



FIG. 2. Diagram yielding the non- analytic momentum 
and frequency dependence of the polarizability no(w,(7) 
in Fermi liquid theory. The solid lines are fermion prop- 
agators. 



Here and below we are use Matsubara frequencies lo = 
2iTnT so that n(w,q) is purely real. n("™f)(w,q) de- 
pends only on the combination fc|| + jk'^ because for 
a circular Fermi line H'^'**"^' (w, q) is a function only of 

|q| -2pF. 

For electrons on a lattice there are additional images of 
coming from fermion transitions with mo- 
menta shifted by reciprocal lattice vectors. The most im- 
portant of these is the transition with momentum G — q 
which might be also close to Q in a typical situation. 
Generally, we expect that the singular contribution to 
fermion polarizability comes from the transitions with 
momenta transfer q and G — q: 

Uo{uj, q) ^ n^'*"^) (^^ q) + nis'^ng) G-q) 

+ n('^™')(tj,q) (9) 
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Here n(°""'^(a;,q) is analytical contribution to fermion 
polarizability coming from fermion momenta far from 
Q/2, its dependence on the momenta and frequency is 
negligible relative to the strong dependence coming from 
the singular parts. 

There are two regimes of q at T = 0. For k 3> AG, 
TIo{lu, q) has a symmetric square root peak at a wavevec- 
tor indistinguishable from G/2. At smaller scales, k <C 
AG, the peaks separate and each peak acquires asym- 
metric form: 11(0, k) ^ -~\fk at fc > and 11(0, fc) ^ 
-|fc|/(AG)i/2 if fc < 0. The qualitative form of n(0,A:) 
is shown in Fig. 3. 



2 ^ ""F 

VP3pFn(0,0) 



For g = (7c, ^(^^jQ) diverges as q ^ Q and w ^ 0. 
For g < Qc the divergence is cut off. It is convenient to 
define a bare dimensionless cutoff, Aq by 



i^o(0,Q)--^ (11) 




-1 -0.5 0.5 1 1.5 



Within the RPA, Aq oc - 9- 

Corrections to the RPA involve diagrams in which elec- 
trons interact by exchange of spin fluctuations. These 
corrections may be organized in a 1/A'^ expansion be- 
cause the spin fluctuation propagator is proportional to 
1/A^ and each fermion loop contains a factor of N . The 
leading diagrams in 1/A^ are the self energy correction 
shown in Fig. 4 and the vertex correction shown in Fig. 
5. These diagrams are infrared divergent; the divergence 
is logarithmic if the external momentum is larger than 
AG and power law if the external momentum is less than 
AG. These two cases require separate discussions. 



III. SCALING AT LARGE MOMENTA 



At large momentum the small difference between the 
Q and G/2 is unimportant and we may write 



FIG. 3. The qualitative form of fermion polarization 
n(0, fc) in RPA approximation, fc = corresponds to the 
commensurate vector G/2, peaks occur at fc correspond- 
ing to q = Q and q = G — q. The inset shows the blown 
up structure of the double peak maximum. 



D{uj,\^) = 



2ttz 



-|- izu: - 



(12) 



After RPA corrections are included, the spin fluctua- 
tion propagator entering the action ^ becomes 

Jo(c.,q)= , \ (10) 

-no(w,q) 

Within the RPA, a second order phase transition occurs 
when the interaction constant 50 controlling the value of 
ZJq ^ is increased through the critical value 



for the spin fluctuation propagator. We use this and the 
fermion Green function (^ to calculate leading correc- 
tions to the fermion self energy, the fermion-spin fluc- 
tuation vertex and the spin fluctuation propagator. We 
find that the fermion self energy and fermion-spin fluctu- 
ation vertex are logarithmically divergent, while the po- 
larization bubble itself which controls spin fluctuations is 
not divergent. We argue that these logarithms sum to a 
power law and we calculate this power law to order 1/7V. 
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FIG. 4. Leading contribution to the self energy. Solid 
line denotes fermion, wavy line denotes spin fluctua- 
tions. 



We begin with the self energy. The leading contribu- 
tion is shown in Fig. 4 and corresponds to 

S(e, p) - .9^ J G{e + lo, V>')D{uj, p - p')dLod?p' (13) 

The form of the fermion Green function implies that 
the energy of the scattered electron is small (~ e) , so the 
frequency transferred to the spin fluctuations is small and 
the scattered electron remains near the Fermi line. Al- 
though p and p' must both be near the Fermi line, the 
angle between them may be large; thus, it will be conve- 
nient to parametrize the momenta p and p' by their polar 
coordinates p, 9, p' , 9' . Moreover, since the momentum 
p — p' transferred to the spin fluctuation is large, we may 
neglect the frequency dependence of the spin propagator 
except as a lower cut off and estimate the large momen- 
tum p — p' neglecting the small differences \p\ — 1 and 
\p'\ — 1. For electrons near the Fermi line the difference 
p — p' is always smaller than 2, so the main contribution 
to the spin fluctuation propagator (|l0| ) comes from the 
term n('*™9)(^^ G - q) in (|), yielding 



Ny/39^ + 36i'2 + 26*61' 



(14) 



The l/\9'\ dependence of D{uj,p — p') at large 9' leads 
to a logarithmic contribution to the self energy and justi- 
fies our assumption that angular deviations are typically 
large. The logarithm is cut off by A which is the largest 
of \e\,p- 1, 9^. Substituting Eq. (|l|) into (|l|) we get 



-IZqC 



ln(l/A) 



(15) 



Here n = 3 is the number of spin components and zq 
is the usual Fermi liquid wave-function renormalization. 
We emphasize again that this expression is only correct 
if the infrared cutoff A is larger than AG. 

Note that S is a function only of energy and 9, so the 
Fermi velocity is not renormalized and the structure (|^) 



that we assumed for the Green function is not changed. 
Note also that the relative value of the renormalization 
I]{e)/{ze) does not depend on the value of the coupling 
constant g or indeed on any other parameter of the the- 
ory except for the number of electron components N and 
spin components n. Thus, even though g is renormalized 
by the interaction, this renormalization is not important 
for the calculation of S(e) and we may expect that the 
logarithms sum up to a power law and that the exponent 
depends only on n and N. We express the renormaliza- 
tion of the self energy as a scale dependent wave-function 
renormalization z{A) and find 



(1/A)" 



(16) 



with a = n /{NVStt) in the large N limit. Note that even 
atiV = 2,n = 3Q;w0.27isa small number, suggesting 
that the leading logarithm approximation is reasonably 
accurate even in this case. 




FIG. 5. Leading contribution to the vertex renormal- 
ization. Solid lines denote fermions, wavy lines denote 
spin fluctuations. 



We now consider the renormalization of the interac- 
tion vertex g. At leading order in 1/A^ this is given by 
the diagram shown in Fig. 5. The evaluation proceeds 
differently from the evaluation of the self energy because 
the momenta of the particle and hole are on the opposite 
sides of the Fermi line and frequency of the spin fluctua- 
tion Lo is not small. It is convenient to use the cartesian 
coordinates p|| andp±, introduced in (0) and Fig. 1. The 
expression corresponding to Fig. 5 is 



Sg 27r(2 - n) f dedp\\dp±_ 



N 



1 



1 



(2^)3 
1 



-P\\ 



Re 



-P\\ 



-P\\ 



This integral is logarithmic. The coefficient of the log- 
arithm may be obtained by scaling e and p\\ by p\ and 
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evaluating the integral over the rescaled e and p\\ numer- 
ically. We find 
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Ntt 



■HI/ A) 



(17) 



with a w 0.75. As in the case of the self energy this 
expression exponentiates, leading to 



5(A) = (1/A)^ 



(18) 



with /3 — a^^^, which in the physically relevant case 
n = 3, iV = 2 becomes (3 ~ 0.08 so the corrections to the 
vertex are very small and we may assume that the one 
loop approximation of the vertex corrections is reason- 
ably accurate in the physical situation. 

Note that both exponents a and (3 depend only on n 
and N but not on any other parameter. This can be also 
seen directly from the action (^ because one can always 
scale away the interaction constant g and the wave func- 
tion renormalization z changing the scales of the S'-fields 
and frequencies. This shows that the effective charge 
controlling the RG flow depends only on the parameters 
n and N and is not renormalized. 

We finally consider the renormalization of the spin fluc- 
tuation propagator D(u!,q). To order 1/7V there are two 
diagrams; a self energy and a vertex correction. Power 
counting shows that the frequency and momentum de- 
pendent terms in Il(uj,q) acquire no additional renor- 
malization beyond the one imposed by the momentum 
dependent z and g, so D{uj,q) is given by Eq. ( [l^ ) 
with z and g replaced by their running values. How- 
ever, there are logarithmic contributions to the mass Aq. 
These come from the diagrams shown in Fig. 6. 



The analytic expression corresponding to the diagram 
in Fig 6a is 



G(6,p + Q)G(e,p)^ 



X G(£-^,p-q)J^(c^,q) ^^J^g (19) 



The self energy part of this diagram (second line in (|1£ 
has an ultraviolet divergence which leads to a trivial shift 
of the fermion chemical potential which we subtract and 
a logarithmic divergence which we obtain by performing 
the integrations in the following order. We integrate first 
over p|| , then over e and finally over p_L, obtaining an in- 
tegral over a;, k\_ and We then find that k\_ can be 
scaled out of the integrals, leading to 



A 



TV 



ca ln(l/A) 



(20) 



Here the coefficient 



dxdib 



X Im 



0.20 



C~^sgn(t^) 
2iL0 ~ 2x + Csgna; 



+ X + iLU + J — X + iuj 



(21) 



where C = v4iw — 4x— T. 

The contribution of the diagram shown in Fig. 6b may 
be evaluated similarly, but in this case no subtraction is 
necessary. We find ultimately a logarithmic divergence 
from the last q± integral. As was the case for the dia- 
gram shown in Fig. 6a, the coefficient of the logarithm 
is independent of g and z: 



SA 



N 



-Cfcln(l/A) 



(22) 





with Cf, w 0.45. By combining the results of these two 
diagrams and integrating the resulting scaling equation 
we find 



A(A) = AoA'' 



(23) 



with exponent rj = 2 ''°""'=^^^" + 0{1/N^). Evaluating 
this formula at = 2 gives f] w 0.15. The formula for 
A (A) implies that the mass term becomes important at 
a scale set by equating the renormalized kinetic term to 
the renormalized mass term, i.e. at 



FIG. 6. Leading contribution to the mass. A, renor- 
malization. Solid line denotes fermion, wavy lines denote 
spin fluctuations. Heavy dot denotes mass operator A. 
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(A)-_yi(A)A = AoA" 



z(A) 

Solving for the scale A we get 



(24) 
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Aa - A = 



(25) 



This result is meaningful only if Aa > AG. At energy 
smaller than A the infrared divergences are absent and 
the renormalization flow stops. One expects that Aq is 
linear in some external control parameter such as pres- 
sure (which would vary the interaction constant) ; thus if 
p — Pc > AG, we would expect that x(0, Q) would vary 
with pressure as 



X(0, Q) cx (p - pc 



2v 



(26) 



IV. SCALING AT SMALL MOMENTA 



angle coordinates used in the discussion of the electron 
self-energy in the previous section. It further develops 
that the dominant contribution comes from processes in 
which the angle 6 pertaining to initial momentum p is 
small relative to the angle 6' pertaining to p'. In this 
limit we may approximate 



2tt 



1 



(28) 



Using this expression for D{lj,c\) and evaluating the di- 
agram shown in Fig 6a gives 



(5A 



2n lln(l/A) 



27TTNh Ai/2 



(29) 



For small energy and momentum the splitting, AG, 
between the peaks in 11(0, q) becomes important. To 
treat this regime it is convenient to expand 11(0, q) for 
momenta small compared to AG. It is also convenient 
to measure parallel momenta in the units of AG and fre- 
quencies in the units of z(AG)AG. In these units the 
fluctuation propagator, 



Z?(co,q) 



27r 



N 



Re\ k 



-11 



h(k 



(27) 



contains a new parameter, 6, which, as we show below, 
controls the renormalization group flow and which itself 
is renormalized. In our units the initial value of 6 is 1/2. 

We analyze the model in the same way as in the pre- 
vious section. The spin fluctuation propagator is more 

singular in the region of large and negative k\\ + when 
the real part of square root in Eq. (|27| ) is small. The 
more singular propagator leads to infrared divergences 
which are stronger than logarithmic, and, we shall show, 
to a first order transition. 

As before, we may consider the renormalization of the 
electron self-energy, the interaction constant and the po- 
larization bubble. Also as before, the renormalization 
of the self energy does not affect the renormalization of 
other quantities; we do not discuss it further. Unlike the 
situation at large momenta, there is no renormalization 
of the interaction constant in the leading order in 1 /N 
because in the diagram of Fig. 5 it is not possible to put 
all fermion lines on the Fermi line and simultaneously 
have the wavy line carry momentum close to Q. The 
leading corrections to both A and b thus come only from 
the self energy insertion in the polarization bubble, as 
shown in Fig. 6a. 

The dominant scattering processes contributing to the 
electron self-energy are those in which the electron mo- 
mentum remains close to the Fermi line. Note that 
for such processes the momentum transfer is such that 

k\\ + ^ ~ 0, so the spin propagator is large. To dis- 
cuss these processes it is convenient to use the radial and 




FIG. 7. Leading contribution to the renormalization 
of the spin fluctuation propagator. Solid lines denote 
fermions, wavy lines denote spin fluctuations. 



The leading correction to b comes from the diagram 
shown in Fig. 7 and may be evaluated similarly. We 
obtain 



5b 
T 



35/6^ 



1 



107V(26)2/3 A5/6 



(30) 



Equation ( p9| ) implies that A decreases exponentially as 
A is decreased. This means that it is not possible to find 
a self consistent solution along the lines of Eq. (^5|) for 
A less than a number of order 1. ■— ■ 

Arguments originally developed by Brazovskiita in a 
slightly different context show that the minimal value of 
A implies a first order transition. The physical reason 
is that fluctuations lead to such a large increase in the 
energy of the critical state that at some point it is favor- 
able to discontinuously open a gap, gaining condensation 
energy and suppressing fluctuations. In the present prob- 
lem, the fluctuations are so strong that for physical values 
of n, N and 6, the flrst order transition happens almost 
immediately as the scale AG is reached and we expect 
that the discontinuities in physical quantities are on the 
scale set by AG. 
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V. PHYSICAL CONSEQUENCES 

In this Section we describe the imphcations of our re- 
sults for observables. We focus on the intermediate scal- 
ing regime discussed in Section III. The observables are 
determined by the spin susceptibility which is restoring 
units: 



52(A) [n(a;,q)+n(^,Q-q) 



^ 27rz{A)vF \ 



\q\ — 2pp iz{A)uj 



PF 



vfPf 



(31) 
(32) 



Here .g^(A) and z(K) are slow power law functions of mo- 
mentum and energy; explicit formulas are given in Eqs. 
(|l^,^8|) . We emphasize that these formulae only apply at 
scales larger than the peak splitting AG. This form for 
x(ti;,q) is essentially the RPA form with small modifica- 
tions due to the momentum and frequency dependence of 

and z. The Eq. ( |3l| ) is written in Matsubara frequen- 
cies. The imaginary part of the analytic continuation of 
x(u;,q) is measurable in neutron scattering experiments. 
The actual form of this imaginary part is somewhat com- 
plicated, due to the structure associated with the bound- 
aries of the particle-hole continuum, so we do not write 
it here. We do discuss in more detail the predictions for 
NMR relaxation rates l/Ti and I/T2, which involve the 
low frequency limits of the real and imaginary parts of 
x(w,q) respectively. 

The static limit of the real part of n(cj, q) is easily 
obtained from Eq. (p2h and is 



j.^ImIIO^ 



AT,vl{\q\-2pFY/^ 



(35) 



We may now calculate the relaxation rates by combin- 
ing the results above with the general relation of the real 
and imaginary part of the susceptibility to the polariza- 
tion operator 



X'(0,q) = 



hm = hm 



<?2[52Ren(0,q) + A] 
g'^lmXl{uj, q) 



--0 Logl[g^ReIi{Q,Ci) + Ay 



and inserting the results into the general expressions for 
relaxation rates 



777— = lim 



jE(Ax'(0,q)f 



where Aq is determined by hyperfine couplings. We see 
that the T2 rate behaves as 



Ren(0, q) 




Ny/POPF 



2'kz{T)vf 



(33) 



There are different regimes for the imaginary part of 
n(u;, q) at small frequency w <C T. We find more conve- 
nient to evaluate /mil directly from the diagram shown 
in Fig. 2 using renormalized Green functions and vertices 
than to analytically continue equation (^2|). The general 
expression. 



lim 



Imn(tj, q) 



/ ImGnie, p + q)ImG7j(e, p) — 
J {2tt 



(Ppdeg'^{e, p) 



(27r)32rcos/i^[e/(2r)] 

is obviously dominated by frequencies of the order of T 
and momenta T/vp so for this calculation we can use 
^^(A) and z(A) evaluated at A = T. If VF{\q\ - 2pf) < 
z{T)T we get 



lim 



Imn(tj,Q) 



Ncpl'^g^T) 
27rzi/2(r)4^^Ti/2 



(34) 



where ci « 0.23. If ?;F(|g| - 2pf) > z{T)T the integral is 
cut off by the external momentum and is 



T2 



(36) 



or, using our previous results a « 0.27, /3 « 0.08 



— =A- BT°°'' 
T2 



Of course, our estimates for a and [3 come from a 1/A^ 
expansion, so for the N = 2 case we may conclude that 
1 /T2 rate is either weakly divergent or non-divergent but 
with anomalously rapid T-dependence at low tempera- 
tures. Evaluating l/TiT similarly we find 



.25 



(37) 



i.e. we expect the l/TiT rate to be weakly divergent 
assuming (as was found at large A^) a > 4/3/3. 
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Finally, we note that this scaling regime only exists for 
energy, temperature and momenta cutoff greater than 
AG. As soon as the relevant scales fall below AG, the 
T — Q transition becomes first order. This transition 
happens at some value, Aq of the controlling parameter 
A. At small but non-zero temperatures the transition 
remains a first order but it occurs at a somewhat dif- 
ferent value of A(r). We may estimate this transition 
line from the argument that if A < Ag, then the or- 
dered phase has lower energy than the disordered phase, 
so /S.E — i?o(A3 — Aq). However, the disordered phase 
has greater entropy, because there is no gap on the Fermi 
line; thus AF = —T^Sq. Equating the two gives 

Tel = To (Ag- A) (41) 

We expect that at scales greater than AG the line of first 
order transitions terminates at a critical point, Tc. We 
estimate Tc{Tc/EfT ~ ^^fAG. 



FIG. 8. Leading contributions to the renormalization 
of the uniform susceptibility. Solid lines denote fermions, 
wavy lines denote spin fluctuations, solid dot denotes 
coupling to the external magnetic field. 



Proximity to the antiferromagnetic transition has also 
an effect pn the uniform susceptibility. We have shown 
elsewheretj that the leading low-T behavior of x(q, 0) in 
the limit q — > is given by 

X[/= lmix(0,k)=^i?(A;,i.)I?(^,q) (38) 

where the coefficient B is given by the diagrams in Fig. 
8, in which fermion mediates interaction between spin 
fluctuations and external magnetic fleld. We denote the 
bare vertex coupling fermions to the external magnetic 
field by ge- 

In order to calculate B we use renormalized 
fermion Green functions and renormalized fermion-spin- 
fiuctuation vertices g. The coupling to the uniform exter- 
nal field is not renormalized; this follows mathematically 
from the fact that in any vertex correction diagram the 
poles in the Green function lines are in the same half- 
plane, so in the limit of vanishing external frequency the 
low-energy contribution vanishes. Alternatively, one can 
show that the absence of logarithmic ge renormalization 
follows from the fact that the logarithmic renormaliza- 
tion given in ( p6| ) does not affect the fermion density of 
states. Estimating the diagrams in Fig. 8 we get 

B ~ z-5/2(e)g2(,),-3/2 (39) 
Combining this with our result for I?(lj, q) gives 

Xu ^ci+ C2T^ (40) 



VI. CONCLUSION 

We have presented arguments suggesting that the 2pF 
density wave transition is first order in c? = 2 spatial di- 
mensions. By way of conclusion we place our results in 
the context of quantum critical phenomena and of the- 
ories of strongly correlated electrons, and discuss some 
implications of our conclusions for the physics of high Tc 
superconductors. 

One may divide quantum critical phenomena in met- 
als into two classes: those in which fluctuations of the 
order parameter are strongly coupled to the particle-hole 
continuum (and in particular may decay into a particle- 
hole pair) and those in which the coupling to fermions 
is irrelevant. If the coupling to fermions is irrelevant, 
then the critical phenomena are described by a theory 
of propagating[-b«srtiis which may be studied by conven- 
tional methodstSllja. If decay into a particle-hole pair is 
possible, then in most cases it is still possible to describe 
the critical phenomena_bjy-^ theory of bosons, albeit with 
overdamped dynamicsBollj. The description in terms of 
a purely bosonic theory is possible in these cases because 
the effect of the critical fluctuations on the fermions is 
small enough so renormalization of the fermions does not 
feed back into the properties of critical fluctuations. 

The one exceptional case is the two dimensional 2pp 
transition. Here the large phase volume available for 
scattering across the Fermi line makes the critical fluc- 
tuations have a strong effect on the electrons; further 
the 2p p singularities in electron response functions mean 
that electrons near the Fermi line which are strongly af- 
fected by critical scattering have a large effect on the 
critical fluctuations. This physics leads to apparently in- 
tractable difficulties in the bosonic model generated by 
formally integrating out the electrons. Specifically, in the 
resulting bosonic model nonlinear terms of all orders have 
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divergent coefRcients and are all relevant in renormaliza- 
tion group sense. Therefore, we used a model in which 
both electrons and spin fluctuations are retained. We 
assumed that the bare propagators and the susceptibili- 
ties have the Fermi liquid form and that the Fermi line 
is not straight near the points connected by the ordering 
wave vector Q. The structure of the 2pF singularities 
enabled us to construct a renormalization group trans- 
formation under which we could treat fermions and spin 
fluctuations on the same footing. 

The solution of the resulting renormalization group 
equations implies that the transition is generically first 
order, because the critical fluctuations are so strong that 
they completely suppress the second order transition. 
There is one special case where the transition is not first 
order. If twice the ordering vector Q is commensurate 
with a reciprocal lattice vector, G, i.e. 2Q — G, then 
the spin fluctuation propagator is less singular, the fluc- 
tuations are weaker and the transition turns out to be 
second order and characterized by the exponents which 
we calculate in a expansion. If |2Q — G| is small, 
the T = transition is ultimately first order but a broad 
scaling regime exists. 

Our explicit calculations were performed for a model of 
a spin density wave transition. Most of our results carry 
over to the charge density wave case, if the number of 
spin components, n, is set to 1. There is one important 
caveat. If the system is sufficiently symmetric, e.g. if the 
Fermi line is circular, cubic terms (forbidden in the SDW 
case by time reversal) may exist in Landau free energy. 
Cubic terms lead to the first order transition but their 
presence would also complicate the analysis given above. 

We now place our results in the context of theories of 
interacting electrons in two spatial dimensions. To pro- 
duce a quantum phase transition one must increase an 
interaction parameter to a critical value which is gen- 
erally large. One must therefore consider the effect of 
the interaction on the non-critical state of the electrons. 
There are two possibilities: one is Fermi liquid theory in 
which it is assumed that perturbation theory may be re- 
summed to all orders. Within this assumption one may 
show that away from any critical point the low energy 
properties are not qualitatively changed from those of 
free electrons, so the transition takes place against Fermi 
liquid background and one may use the theory given here 
to calculate the extra effects due to criticality. 

An alternative possibility is that Fermi liquid theory 
is not a correct description of the low energy properties 
even far from criticality. For example, it has been ar- 
gued that in strongly correlated two dimensional^wiels 
the spin degrees of freedom form a "spin liquid"Ot3. A 
spin liquid possesses a Fermi surface, spin 1/2 fermionic 
excitations with constant density of states at low ener- 
gies and a particle-hole continuum but the fermions inter- 
act via a singular interaction which, anmng other things, 
changes substantially 2pp singularitiestj. The singular 
interaction causes anomalous temperature dependence of 
susceptibility and NMR relaxation rates even away from 



the criticalitytS and changes the critical properties. The 
spin-liquid to-antifcrromagnet transition has been stud- 
ied elsewherell3. 

The high temperature superconductors are quasi two 
dimensional materials with spin dynamics which have 
been claimed to be controlled by a quantum critical fixed 
point. It is still controversial whether Fermi liquid the- 
ory is the correct starting point for a description of the 
low energy electron physics. It is therefore interesting to 
compare our theoretical results to the known magnetic 
properties of high T,. materials. 

For a two dimensional Fermi liquid the bare polar- 
izability is so strongly peaked at 2pp that it is most 
natural to assume that the transition occurs at Q = 
2pp. This is consistent with neutron scattering data on 
La2-xSrxCuOi^ in which strong and temperature depen- 
dent peaks were observed. These peaks are centered at x- 
dependent wavevectors Q(a:) which are claimed to be 2pp 
wavevectors of the LDA band structureE£l. Therefore if a 
Fermi liquid picture is appropriate, it is natural to expect 
our results to describe experiments on La2-xSrxCu04. 

In fact, they do not. There is no sign of a first or- 
der transition; the magnetic properties apparently evolve 
smoothly with doping. It is conceivable that disorder due 
to random positions of the Sr dopants masks the first 
order transition. However, for all Sr concentrations in- 
cluding X ~ 0.14 there is a wide temperature regime in 
which the NMR relaxation rates and uniform susceptibil- 
ity vary with temperature. Roughly, the coppei:JL/(riT) 
and I/T2 are— inversely proportional to l/7oEJ while 
Xu ^ A + BTE3. These temperature dependences are not 
consistent with our results, Eqs (|^,||,|o|). We conclude 
that the magnetic properties of La2-xSrxCuOi are not 
well described within a Fermi liquid approach. Two al- 
ternatives have been proposed: one is that tiie fermions 
are in "spin-liquid" regime described aboval3; another 
is that the critical behavior is due to prottagating spin 
waves only weakly coupled to the electronaJ'O. 

Another class of materials to which our results might 
be relevant are the low dimensional organics such as 
TTF-TCNQ. These materials have strongly anisotropic 
transfer integrals ia 3> it S> ic, leading typically to an 
open Fermi line, but with non-negligible curvature. The 
curvature implies that the physics of these materials is 
not strictly one-dimensional and makes it possible that 
the theory developed here is relevant and explains the ex- 
perimental observation that the spin density wave tran- 
sition prOjduced by lowering the temperature is weakly 
first ordeJi^l. 

Note added. After this paper was completed we re- 
ceived a preprint from A. V. Chubukov analysing the spin 
density wave transition with Q = G/2 but Q ^ 2pp. He 
found a logarithmic renormalization very similar to the 
one we found in the intermediate scaling regime discussed 
in section HI and showed that this implies that the ex- 
ponents characterizing the critical point he analysed dif- 
fer from the exponents characterizing the general case 
in which Q 7^ G/2 and Q, ^ 2pp analysed by previous 
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